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Abstract. We show how the correction to the calculation of the mass in the original 
relativistic model of a rotating star by Hartle [B], found recently in m, appears in the 
Newtonian limit, and that the correcting term is indeed present, albeit hidden, in the 
original Newtonian approach by Chandrasekhar |2j. 


1. Introduction 

The original treatment aimed at the study of (rigidly) rotating stars in a perturbative 
scheme is due to Chandrasekhar for Newtonian gravity, back in 1933 [2]. It was only in 
1967 that Hartle put forward the model within the realm of General Relativity [6]. 
Although the study in [6] covers any barotropic equation of state, the work in |2] 
focuses, from some point onwards, only on polytropic equations of state, i.e. of the 
form p = Kp 1+ ~ for some constants K and n, where p and p denote the pressure and 
the mass density of the star. The relationship between the Newtonian and the GR 
approaches was presented in [6], and the GR procedure was found to be consistent with 
the Newtonian case by taking care of the suitable limit. 

However, it has been found recently [10] that the computation of the total mass 
of the rotating configuration as a function of the central density shown in [6] has to be 
amended by a term proportional to the value of the background energy density at the 
surface of the star. That value is zero for certain equations of state (including polytropic 
EOS), but it does not vanish necessarily (for instance in models of strange quark stars 
0)- At the time of writing [TO] we did not pursue the study of the Newtonian limit, since 
no such term seemed to appear in the literature concerning Newtonian stars. However, 
as we show next, that term contributes to the Newtonian limit, and appears indeed, 
although implicitly, in the original work by Chandrasekhar [2]. appearance of that term 
had been somehow forgotten, even in the review of the Newtonian approach in [6]. 
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We first concentrate on the computation of the mass as stated in [2], expand that for a 
general case, and show how the expressions in PJ for polytropic equations of state follow 
indeed. In page 396 of [2] the mass is claimed to be given by 


M = 2n 


j J pr 2 drdp, 


( 1 ) 


where p := cos 6, and {r, 6, are spherical coordinates, so that r and 6 are the radial 
coordinate and azimutal angle on the sphere respectively. 

In agreement with the next equation in [2], as we will show later, © stands for 
the integral over the deformed volume. Indeed, the shape of the star is described in [2] 
to be the sphere of the background configuration plus a deformation at first order in a 
perturbation parameter v, which corresponds to a second order in the angular velocity 
( jj 2 over the value of the central density, see (10) in [2]. The value of the central density 
is denoted by A in [2], but we will use p c here. 

Let us first review just the necessary of the Newtonian treatment in order to obtain 
the expression of the total mass of the rotating star, suitable to be computed by solving 
the relevant problems at different orders in the perturbation. We follow essentially the 
description of the Newtonian approach as made in ( 6 ], and will compare with that in [2] 
when necessary. 

The non-rotating (spherically symmetric) configuration is described by the mass 
density p(°)(r), pressure p(°\r) and Newtonian potential U^(r). The radial variable r 
runs from 0 to a in the interior and r > a corresponds to the vacuum exterior. The three 
equations of structure that govern the configuration are a barotropic equation of state 
pO) — p(°) (p(°)), the hydrostatic equilibrium first integral and the Poisson equation, i.e. 

r l dp^ (s' 


7 = 


-ds + f/ ( 0 ) (r), V 2 C ( 0 ) (r) = 4 tt Gp ( 0 ) (r), 


( 2 ) 


/o p(°)(s) ds 

where the constant 7 is identified as the chemical potential. We use V 2 for the flat 
Laplacian in spherical coordinates {r, 9, f>}. In spherically symmetric configurations the 
mass function is given by 

M(°)(r) = inf p(°\s)s 2 ds. (3) 

Jo 

Given that M^(0) = 0 and that regularity at the origin r = 0 implies dU^ / dr ( 0) = 0, 
the mass and the potential are related by 

dU^Xr) G A/r(0) ,, 

- -±± = —M^\r). 

dr r z 

The system of three equations can be integrated in terms of boundary conditions at the 
origin and thus provide, e.g. , the total mass of the star, Mg := M^°^(a), in terms of 
the central density p c = p^( 0). We denote that function by Mg ) \p c ). 

Consider now the (perturbed) rotating configuration. The mass density of the 
rotating configuration p(r, p) is expanded perturbatively to first order in v as 

p(r,p) = p ( 0 ) (r) +np ( 2 ) (r, p) + 0[y 2 ). 
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A new radial coordinate R is now chosen so that it labels surfaces of constant density 
in the rotating configuration by (see 0) 

The interior of the rotating star is therefore defined by R G (0, a ) by construction, and 
its surface located at R = a. The change between R and r must thus have the form 

r(R,p) = R + v((R, /i) + 0(v 2 ) (5) 


for some function f(R,p), which thus describes the deformation of the surface [2, [ 6 ] . 
Before considering the equations governing the rotating configuration, in the next 
section, let us expand the integral ([I]), which reads explicitly 

r 1 ra+v£(a,/i) , . 

(p(°)( r ) _|_ V pG)(j., pyj r 2 drdp + 0(v 2 ). ( 6 ) 


M — 2n 


-1 jO 


In [2] the radial coordinate r is, instead, conveniently rescaled to a new radial coordinate 
f (eq. (9) in [2]). The shape of the star is then described in |2] to be the sphere of the 
background configuration £i plus a deformation denoted by df\ at first order in v. For 
a polytropic equation of state P = Kp 1+ ™, is the first zero of Emden’s function with 
index n, denoted by 9(f) in [ 2 ], and the deformation corresponds to df(f i,/i) for some 
function df(f, /i) that can be extracted from the terms in the v factors in equations (36) 
and (38) in [2j. Obviously v((R, p) scales to df(f,p) as r scales to f (9) in |2]. Note 
that df contains v. After using (7), (9), (36) in |2] for the polytropic equation of state, 
m can be shown (see below) to translate, up to order v, to 


M = An 


n 


AnG 


A n 


3/2 


A 


Ri+d£i 


(9 n + vne n -^ 0 )f 2 df, 


(7) 


as it stands in page 396 in J2J, where dfi denotes the l = 0 part of df(f i,/i), which 
equals dfi = — u , 0 o (£ 1 )/6 ), (£ 1 jj], see (38) in [2], where f>o(f) is an auxiliary function that 
satisfies (37i) in |2J. df i is the expansion of the star, as noted in [2]. Only the / = 0 
sector contributes to the integral. 

In order to obtain (JTJ) and go further let us develop ([ 6 j). Since the Jacobian of the 
change (JSJ) is 1 + vd( / dR , the integral ([ 6 ]) expands as 


M = 2n 


I I p(r(R, /i),/i)(A 2 + 2v(R)(l + v^-)dRdp + 0(v 2 
J -i J o oR 


= 2n 


= 2n 


= 2n 


' rp {0] (R)(R 2 + 2vfR)(l + v^OdRdp + Q(v 2 ) 
-i Jo oR 

d( 


rj; 

u: 


p w (R)R 2 + v 2p {0) (R)(R + R 2 p {0) (R)^r 


pW(R)R^ + v (-R\^- + 


dpd 


dR / 
(r 2 p^\r)C) 


An j a p {0) (R)R 2 dR-2nv f 
+0(v 2 ) : 


-1 JO 


dR OR 

R 2 ( P dRdp + 2nv 
dR 


dRdp + 0(v 2 ) 
dRdp + 0(v 2 ) 
a 2 p {0) (a)((a,p)dp 


( 8 ) 


| Since 0'(£i) < 0 [2], — 0'(£i) always appears as |0'(^i)l i n E- 
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where the relation (J3J) that defines R has been used in the second equality. The first 
term in the final expression (jS]) corresponds to M^\a) by (J3]h The second term is more 
easily recognised by developing the relation (j3J) as follows. Given that for any /(r,/i) 
we have 

f(r(R, /i), /i) = f(R, n) + vf(R, n)C + 0(v 2 ), (9) 

where the prime denotes differentiation with respect to the first (or only) argument, the 
relation (j 3 D provides, in particular, 

dp(°\R) 


p( 2 )(f?,/i) = -C0R,h) ; 


dR 


( 10 ) 


The expression (jSJ) can thus be also written as 


M = 4vr I" p (0) (R)R 2 dR + 2nv [' p {2) (R, p)R 2 dRdp 

Jo J -i Jo 

a 2 p ( 0 ) (a)C(a, p)dp + 0(v 2 ). 


2nv 


( 11 ) 


From now onwards let us denote by a f 0 (subindex 0 ) the part of any function 
/ parallel to the Legendre polynomial P 0 (p)(= !)• I 11 other words, /o(-) : = 
1 //(-, p)Po(p)dp. We will also refer to /o as the l = 0 sector of /. The mass (fill) 
thus reads 

M = An f p ( '°\s)s 2 ds + Anv [ p^\s)s 2 ds + Anva 2 p 1 ' 0 ' 1 (a)( 0 (a) + 0(v 2 ). ( 12 ) 

Jo Jo 

The fact that only the / = 0 sector contributes to the integral is now explicit. 

For polytropic equations of state, after using (7), (9), (36) and (38) in [2], (11211 
directly translates, up to order v, to 


M = An 


(n 


AnG 


A n 3 


3/2 


n 2 ^ + vf' ne*- l wk |. (13) 


fl'Ki 


which is not difficult to show to be equivalent to (J7J) irrespective of the equation that 
the function 6(f) satisfies. 

The crucial point here is, let us recall, that the function 0(r) is Emden’s function, 
for which 6 ) (^' 1 ) = 0 by construction, which is equivalent to p^(a) = 0. The above 
expression (TT3|) for the total mass is obviously presented in [2] without the last term, 
which vanishes (see above (40) in [2]). However, in general, the mass density p^(R) 
of the background spherical configuration does not have to vanish necessarily at the 
boundary R = a. The expression of the total mass in [2], made explicit for a class of 
equations of state for which the mass density vanishes at the surface of the star, seems 
to have misled many authors to forget the third term in dT2]h Even the author himself 
forgot, many years later, to include that term when exploring homogeneous (constant 
p) stars in GR [H] . The correction to the calculation of the mass of homogeneous stars 
can now be found in [9] . 

The third term in (fT2]h proportional to p^°\a), corresponds, precisely, to the 
Newtonian limit of the term found in [TO] that amends the “change in mass” computed 
in [ 6 ]. That is shown in the following section, where we very briefly review the equations 
for the perturbed configuration needed in both Newtonian gravity and GR. 
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3.1. Newtonian gravity 


We only need to consider the l = 0 sector of the perturbation. As in the background 
configuration, apart from the given barotropic equation of state, the perturbation at 
first order in v is governed by a hydrostatic equilibrium first integral and a Poisson 
equation 


0\R) + UR) 


dU<®(R) 2nGp c R 


_[ n2 

R 2 dR 


dR 

dU { »\R ) 
dR 


= -47rGCo(R) 


= 0 , 

dp<®(R) 
dR ' 


(14) 


(15) 


where the Poisson equation for the nonrotating potential has been used in the second 
equality. Note that from (HUD we have p^(R) = —Co (R)dp^(R)/dR, so that the right 

/ n\ 

hand side of (fl5D can be also expressed as 47rGpg (R). It is important to note that 
the domains of definition of these equations are given by R e (0, a) for the interior and 
R > a for vacuum, and suitable boundary conditions (including a regular origin and 
asymptotic flatness) are imposed accordingly. 

It is convenient to change the functions {Uq y Co } that describe the configuration 
to a new set {M^fp^}, suitable to be compared with the relativistic model, defined as 
follows, 


( 2 ) | 


M' 2 >(fi):=47r| o pg 

m) - <*). 


s)s 2 ds = —47r 


rR dpf"Hs ) , 

Co S s ds : 

o ds 


(16) 

(17) 


The definition can be expressed in terms of the pressure and density of the 
background configuration by differentiating the hydrostatic equilibrium first integral 
for the static configuration (first equation in (J2D), which provides 
dUl°\R) 1 dpW(R) 


dR 


pW(R) dR 


= 0 , 


so that 


UR) = -p m (R) p’o( R )- (18) 

On the other hand, the second order Poisson equation (fTU]) can be expressed in terms 
of the pressure perturbation factor by using (fT8]i to get 


1 d ( D2 dG\R) \ _ .(0, 


R 2 dR R 


dR 


(19) 


We can also rewrite the expression for M^ 2 \ (TT6|) . using (TT8|) and f[I7|) . which in 
differential form reads (see (15) in [U]) 


dMUR) = ^ R2 djU, >m 


dR 


dp(°) 


( 20 ) 
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The equation for p* Q is obtained as follows. Combine (fT9l) with (l20l) to get rid of 
Pq and integrate once taking into account that M <k2 \ 0) = 0, by construction, and 
dU^/dR | R=0 = 0 for a regular origin. We thus obtain 


= g (2) 

dR R 2 y h 

in anlogy with the background configuration, 
hydrostatic equilibrium first integral ffT4j) 

iUg\R) + d_ A (R) dVpW \ 
did \ dR J 


( 21 ) 

Finally, take the derivative of the 


AirGpc 


dR 


dR 


R = 0, 


( 22 ) 


and use (12TD and (TTT1) to obtain (see (15) in [5]) 

d jm = G M v ( R )+ ^ R . 

dR K z o 


(23) 


The system of equations for the functions {M^fp^} is formed by (I2U]1 and (ITIli on the 
domain R G (0, a). As in the background configuration system, this problem allows 
us to integrate {M^ 2 \pl} given boundary conditions at the origin. In particular one 
can compute Mi?' 1 := M^ 2 \a) as a function of the (total) central density p c , and thus 
construct a function M$ (p c ). In order to add this function to the contribution from the 
background configuration Mf\pf) it is, of course, necessary to choose p c = p c , so that 
p c becomes a parameter of the whole perturbed configuration. That implies choosing 


Po(0) = 0. 

The total mass of the rotating configuration (1T2T) . taking into account (|2D, (ITBD and 
mi, can be expressed as 


M = M (0) (a) + uM ( 2) (a) + 4vru 




p {0 \a)p* 0 (a) + 0(v 2 ). 


(24) 


Note, again, that this sum makes sense once the functions involved are computed given 
common boundary data, in terms of a common set of parameters, as for instance p c . 
Nevertheless, the choice of parameter used to compute those functions is irrelevant for 
our purposes. The contribution of the perturbation to the total mass in Newtonian 
gravity is given by 

m t ) = M { 2 \a ) + G M(o)^ P {0) ( a )Po( a )- ( 25 ) 

The second term in the above expression is missing in the first equality of equation (18) 
in [6]. 


3.2. General Relativity 

The metric up to second order in some parameter e, following the notation in )6] (see 
also [ID]), reads 

g £ = - e^ r) (l + 2 e 2 h(r, 9))dt 2 + e A(r) ( 1 + 2g 2 — ^ 9 ^ ) dr 2 

\ r — 2 M ) 

+ r 2 (l + 2 e 2 k(r, 9 )) (d9 2 + sin 2 9(dp — eu(r)dt ) 2 ) . 
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Given that we follow |6] and |T0]J in this section, we also use geometrized units for 
convenience, so that G — c — 1 unless otherwise stated. We can fix the (dimensionless) 
perturbation parameter e in analogy with the formalism developed in [2] for the 
Newtonian model. To this aim we set e 2 = v = uj 2 /2nE(fS), where 77(0) is the energy 
density of the background configuration at the origin, and c o is the constant angular 
velocity of the fluid, as in the Newtonian treatment. Therefore, the quantity that drives 
the perturbations in |6j is expressed here by Gt H = ^j2irE(0)v, whereas the constant 

used in na is identified with in this convention^ 

We shall keep the perturbation parameter £ and the constant in this section in 
order to ease the comparison with [10j, although the identifications will be made explicit 
when the Newtonian limit is taken. 

As in the Newtonian case, we only need focusing on the l = 0 sector of the solution 
for our purposes. The coordinate r is fixed by choosing k 0 (r ) = 0 pa] (see also [ID] for 
a discussion on the choice of gauges). The asymptotically flat vacuum solution is given 
by [6] 


^vac{p) 


2 _ _ e ~Kac(r) , ,vac( 


LU 


2 J 
—fN ’ 


5M J 2 

r — 2 M r 3 (r — 2 M) ’ 


m v 0 ac (r ) = 5M 


E 

r 3 ’ 


(26) 


where M, J and 5M are constants. In the analysis of the background and first 
order configurations, M and J are identified as the background mass and the angular 
momentum, respectively. The equations governing the background and first order 
configurations are used to compute M and J given suitable data at the origin. We 
refer to [6j for a full account (see also mini). The constant 5M, still to be determined, 
is identified with the “change in mass” due to the second order perturbation, or simply 
the contribution to the mass at second order, since the angular independent part of g tt 
is given by 

Qttlr^oo = 1 - ~{M + £ 2 6M) + O ^ . (27) 

The l = 0 sector of the (second order) perturbation interior configuration is completely 
determined by the pair of functions {m 0 (r), Pf (r)}, with 


D*_ 

-t n • — 


P( 


( 2 ) 


2 (E + py 


(28) 


where E and P are the energy density and pressure of the static background interior, 
respetively, and P^ (r, 9) the perturbation to the pressure (see nm for this alternative 
definition of the same function in [6]). The system of equations that {mo, -P 0 *} satisfy 
are to be fulfilled in the domain r G (0, a), with suitable boundary conditions, and read 


§ This value is related to the usual choice in numerical works (see e.g. [5]), in which an estimate of the 
mass shedding frequency is chosen, say, ft* = i J M/a 3 . For a constant density star it is easy to check 
that ^2nE(0) = y/3M/2a 3 = . 
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dm o 
dr 

dPl 

dr 


. 2 / t-, dE 1 a (du 

Anr 2 (E + P)-P n H- Wr A — 

V J dP 0 12 J \dr 

(E + P)r 2 r 2 m 0 

71 r-2M 0 (r - 2M) 2 


2 

- 7T r “J 


3 , d'i 


dr 


UJ 


1 r 4 j 2 f du\ 
12 r — 2M y dr / 


Id/ r 3 j 2 u 2 \ 

3 dr yr — 2 M J ’ 


(29) 

(30) 


where u)(r) := u(r) — and j(r) := exp[— {y + A)/2]. 

The value of 5M is determined in terms of interior quantities using the matching 
conditions for the exterior and interior problems to second order provided in |10j. In 
particular, a function m 0 (s) for s G (0, oo) constructed by joining m 0 (s) and m“ c (s) 
across s = a is not continuous in general, since it presents a jump proportional to E(a). 
The result is m 

T 2 a — 2 M 

5M = mo(a) H—- + 47ra 3 ———P(a)P 0 *(a). (31) 

The third term accounts precisely for the jump in mo. Let us stress here that this term 
does not appear in the analysis in [6] because the function m 0 (s) is assumed, a priori , 
to be continuous across s = a. The relationship between the continuity of mo and the 
vanishing of E(a), and therefore the validity of the assumption for certain equations 
of state (in particular the polytropes) is proven and discussed in [TU] . using the fully 
consistent theory of perturbed matching to second order in [8]. 

As in the Newtonian case, the background quantities E(a), M and J, and the 
perturbation ones, m 0 (a) P 0 * (a) are to be computed by solving the corresponding system 
of equations given the (common) relevant data at the origin. In [6] the parameter chosen 
is the central desnsity p c , but, as mentioned above, that choice is not relevant for this 
discussion. 


3.3. Newtonian limit 


Our purpose now is to obtain the Newtonian limit of 5M in (13T| and compare it with the 
contribution to the mass of the perturbation in the Newtonian approach, Mf ; , given by 
(1251) . First, though, it is convenient to find the Newtonian limit for the system (1291) and 
(1301) in order to relate {mo, Pf } with the pair from the Newtonian approach. 

This is achieved by performing an expansion in powers of 1/c as (see |6j) 


M = 


G 


Mf ] + O 


E(r) = 


G 


p(°)(r) 


+ o 


P(r) = 


G 


p(o) ( r 'j 


+ o 


u(r) 


dNER + o (I) 

r \ J 


Po(r) = ^ P * 0 { r ) + O Q) , m 0 (r) = ^ m 0 (r ) + O (J^) , 

for some functions m 0 and pg, where p^ (and p c ), p^ and correspond to the 

functions describing the Newtonian background configuration. Note that, concerning 
the first order (in e), the function u)(r) is constant at lowest order in 1/c [B]. Given the 
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system (f29j) and (1301) . the pair {m 0 ,Po} thus satisfies (see (102)-(103) in [6] ) 


dm 0 

dr 

dp*o 

dr 


= Airr 2 


dp® 

dp(0 

Grhi 


PPm 


o . 4vrGp c 
- H-z— r. 


(32) 

r- 3 < 33) 

Compare this system of equations with (120]) and (l23j) . The functions arising from the 
newtonian limit {mo,Po} an d the functions in the perturbed Newtonian model {M^ 2 \p^\ 
satisfy the same equations in the same domain r, R e (0, a). Therefore, the pair {m 0 ,p^} 
is equivalent to {M^fp^} for r,R < a. We can now substitute {moJol by {M^fp^} 
in the following. 

The Newtonian limit for (13TT) is obtained following procedure above together with 

O ( dl , (34) 


(35) 


_ 

5M = 4 5M 


for some 5M , from where (13TT) becomes 


5M = M (2) (a) + An a ) P {0) (a)p* 0 (a) + O . 


GM 

Comparing this expression with ([25]) we finally find 


5M = Mr 


( 2 ) 


T > 


that is, the Newtonian limit of the contribution (to second order) of the perturbation to 
the mass in GR, found in [TD|, is non zero, and agrees with the same quantity computed 
in Newtonian gravity. 


4. The Newtonian matching conditions 

As a final remark, let us comment on the boundary conditions at the surface of the 
star, the matching between the interior and exterior problems at each order, involved in 
the Newtonian approach. Some objections to the Newtonian matching problem stated 
in [2] were raised in jTj. Those were finally solved by Chandrasekhar and Lebovitz in 
[i.3] by properly formulating the matching and producing the same results. However, 
[3] concerns, again, only polytropic equations of state, and the matching conditions 
are obtained only for that case, which in particular satisfies p(a) = 0. Let us, for 
completeness, deduce the matching conditions for the perturbed Newtonian potential 
in the general case, which, as expected, turns out to be compatible with the obtaining 
of the perturbed mass (124|) . 

Consider any function / that depends on v on two arguments by f(a+v((a , p), /i, v). 
Let us use d v to denote a derivative with respect to the third argument, and a prime 
' with respect to the first, and define /^°^(a,/r) := f(a + v((a,p),p,v )|„ =0 = /(a,yU, 0), 
and /^(a,/i) := d v f(a + v((a, p), p,v)\ v=0 . Assume now that / satisfies the equation 

f(a + v((a,p),p,v) = 0. 
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Evaluating the equation at v = 0 we obtain 

/ ( °W) = 0, (36) 

while differentiating with respect to v, and then evaluating at v = 0 we get 

/ (2) (a,/i) +/ (0) '(a)C(a,/i) = 0. (37) 

Consider now two families of interior and exterior problems, defined by some 
parameter v, with their respective Newtonian potentials Ui nt (r, p,v) and U ext (r, p,v). 
The matching of the problems at each v accounts for the equality of the potentials 
and their radial derivatives at the common boundary, the matching surface. Let us 
define the matching surface by r(a,p) = a + v£(a,p). The conditions are, therefore, 
Uint(a + v((a,p),p,v) = U ext (a + v((a,p),p,v) and U' nt (a + v((a, p), p, v) = U' ext (a + 
v((a, /i), /i, v). Let us finally use the notation [g\ := g int — g ex t, so that the matching 
conditions read [U](a + v((a, /i), /i, v) = 0 and [U'](a + v£(a, p), p, v) = 0. These two 
functions satisfy the requirements for / above, so it is now just a matter of applying 
equations (j36j) and (1571) for both [U] and [[/']. The four equations thus obtained read 

[C/ (0) ](a) = 0, [t/ (o)/ ](a) = 0, 

[U {2) ](a,p) = -[C/W'](a)C(a,/r) = 0, [t/ (2)/ ](a, p) = -[£7 (0) "](a)C(a, p), (38) 

where we have used that the background potentials t/^’s do not depend on /i. The last 
equation in (1581) yields 

[U {2) '](a,n) = -47rGp (0) (a)C(a,/i). (39) 

As expected, although the background potential, its first derivative and the 
potential at first order are ‘continuous’, the radial derivative of the Newtonial potential 
at first order in v suffers a jump, proportional to The perturbed mass can now 

be computed from the Newtonian potential, and it is straightforward to show that this 
jump generates the term proportional to p(°\a) in (15T)) . 

Finally, it can be shown that the matching condition for h! in [HI], which suffers a 
jump proportional to the jump of m, agrees with (158j) after taking the Newtonian limit. 
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